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1 Introduction 


The Acoustical Theory for Design of Aircraft Cabins (ATDAC) is a computer pro- 
gram developed to predict interior noise levels inside aircraft and to evaluate the effects 
of different aircraft configurations on the aircraft acoustical environment. The primary 
motivation for development of this program is the special interior noise problems asso- 
ciated with advanced turbo-prop (ATP) aircraft where there is a tonal, low frequency 
noise problem. This report describes the theory and equations implemented in the AT- 
DAC program. A description of the program structure and the operating instructions 
is contained in the user’s manual: “Interior Noise Prediction Methodology-ATDAC 
User’s Manual.” [l] 

Prediction of interior noise levels requires knowledge of the energy sources, the 
transmission paths, e.g., structural elements and the relationship between the energy 
variable and the sound pressure level. The energy sources include engine noise, both 
airborne and structure-borne; turbulent boundary layer noise; and interior noise sources 
such as air conditioner noise and auxiliary power unit noise. Since propeller and engine 
noise prediction programs are widely available, they are not included in ATDAC. Air- 
borne engine noise from any prediction or measurement may be input to this program. 
ATDAC includes modules which predict the contribution of boundary layer noise and 
structure-borne noise from the engine. Several ways of predicting structure-borne noise 
are discussed in Section 4. The turbulent boundary layer noise algorithm is presented 
in Section 5. The acoustical power from the noise sources inside the aircraft may also 
be input to ATDAC. 

The energy produced by exterior noise sources propagates to and impinges on a 
panel at some point. The transmission loss module predicts the energy behavior at 
the panel by predicting the transmission loss and reflection co-efficients of the panel. 
The theory to predict transmission and radiation properties is described in Section 3. 
The transmission loss and reflection co-efficients are put in a database of panel designs. 
Any experimental panel data may be included in the database. 

The energy sources and panel properties are combined with a description of the 
aircraft s geometrical configuration to predict the interior noise levels. A geometry 
module defines the volumes and panels for a particular aircraft. The power balance 
takes the source data, the transmission and reflection data and the geometry data and 
solves for the net intensities into each volume. Then the interior sound pressure levels 
are calculated. The form of the power balance is based on the room equation and is 
described in the next section. 
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2 Energy Balance and the Room Equation 


2.1 Introduction 

The ATDAC noise prediction method is based on a modified form of the acoustic 
room equation. The acoustic room equation, in its common form, was published by 
Embleton [2], The room equation has been widely used, especially in the field of 
architectural acoustics. Many empirical corrections have been proposed for the room j 

equation but none will be used in this implementation so that the method will be as 
general purpose as possible. The room equation, however, has been modified for use 
in the present noise prediction methodology. The choice of using the room equation I 

represents a compromise between accuracy and simplicity. 1 

The accurate prediction of sound contained in a multiply- connected system of en- I 

closures is a difficult problem involving the knowledge of both the structural and the \ 

acoustic behavior and how they interact. For a system of any size at any frequency ! 

range but the lowest, the set of equations resulting from numerically solving the wave { 

equation in the cavity coupled with the equations of motions for the structure is very 
large. When the size of the system of equations is combined with the extensive mod- 
elling that must be done to generate those equations, it is prohibitively expensive to 
analyze most realistic structural/acoustic systems utilizing first principles. This is espe- 
cially true when multiple trade studies must be done to evaluate many alternatives. By I 

using the room equation, the sound pressure levels in large structural/acoustic systems ! 

can be evaluated quickly. Because the room equation uses an energy or power balance 1 

formulation, the structural equations are decoupled from the acoustical equations and 
therefore detailed analysis may be done on specific portions of the model and later the 
results can be substituted into the overall prediction. In the same way, experimental 3 

results for specific components may be used. j 

2.2 Assumptions 

In order to use any theory well, the assumptions and limitations must be understood. 

This is particularly important in this case because the simplicity of its use may lead to 
trying to draw more conclusions than are appropriate. 

These assumptions are as follows: ! 

• the system has reached steady state conditions. 1 

• the sound field in the room is completely diffuse and reverberant so that the sound : 

energy density is uniform throughout the room. This implies that the absorption ' 
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is small and uniformly distributed. 


• the mean free path (MFP), which is the average distance sound travels between 
reflections, must follow 

4V 

MFP = — , (1) 

where V is the room volume and S is the total surface area of all the walls. This 
criteria is usually met if the lengths of the sides of the room are all of the same 
order of magnitude. 

• All noise sources are sufficiently distributed so that the direct field may be ignored. 

Because the major sources are spatially distributed in aircraft, the assumption of 
a diffuse field is reasonable except where the response is largely dominated by one or 
two modes. 


2.3 Derivation 

In terms of engineering units, the room equation can be written 

= + < 2 > 

where p is the acoustic pressure, W is the power generated in the room, p is the air 
density, c is the speed of sound in air, Q$ is the directivity factor of the source, r is the 
distance from the source, and R is the room constant. The term represents the 
energy in the field directly radiated by the source and the term ^ represents the energy 
in the reverberant field. Splitting up the energy terms into direct and indirect fields 
requires detailed knowledge of the energy sources. A modified room constant can be 
defined in which the term describing the direct field is absorbed into the reverberant 
energy term giving the equation 


W* = 


(3) 


where R' is the modified room constant. After absorbing the direct field term into the 
reverberant field term and realizing that the intensity in a completely reverberant field 
is 


jpl 

4 pc 


the room equation can be written as 


R'l = W. 


(5) 
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Because the room equation was derived for a source in a room, the input power is 
modified in this derivation to include all the power coming into a room. By expanding 
the definition of W this way, equation (5) can be rewritten 

#<?/ Surf aces #of Sources in Room 

R'i = £ w t + £ w h ( 6 ) 

*=l 1=1 

where W k is the power coming into the room through the kth surface and Wi is the 
power generated by the Ith source in the room. The power coming through a surface 
is given by 

w k = 1ST, ( 7 ) 

where I is the intensity incident on the other side of the surface, S is the area of the 
surface and t is the transmission co-efficient of the surface. If the room is connected to 
other rooms, some with known intensities (type A) and some with unknown intensities 
(type B), the room equation may be rewritten as 

if of Type B Surfaces if of Type A Surfaces if of Sources in Room 

R I — I k SkTk — I m S m T m + Wi. (8) 

fc=l m= 1 l = i 

On the left hand side of equation 8 all the intensities, /, are unknown but on the 
right hand side all the variables are known. As it is, equation 8 is ill-posed but more 
equations can be written. If one equation is written for each room or volume where the 
intensities are not known and the intensities are known outside of all exterior walls or 
surfaces, then a system of equations can be written in matrix form for a set of n rooms 
or volumes 


(SYM.) 


where 


it-of Exterior Surfaces # 0 / Sources in Volume i 

Wi = IimSimTim + ]T W it 

m=l ;=i 


Srij - 


#of Sur faces connecting Volumes i&j 

RkT k . 

k = 1 
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Every term in equation (9) has been defined except R', the modified room constant. 
The room constant is typically defined as [2] 


R = 


Sta 

1 — a 


( 10 ) 


where St is the total surface area of the boundary of the volume or room and a is the 
average statistical absorption co-efficient. The 1 — a term represents the ratio of energy 
left after the first reflection or the ratio of the energy supplied to the reverberant field 
from the direct field. In order to make the room constant account for ail of the energy 
lost from the room, the energy lost during the first reflection must be included. The 
modified room constant which includes all the energy is 

R' = S t a. (11) 


The average statistical absorption co- efficient is defined by 


a — 


Boundary Surfaces 
Boundary Surfaces q 


( 12 ) 


The statistical absorption co-efficient a is defined as the ratio of the power absorbed 
(n„) to the power incident (Ili) 

(13) 


a — 


n. 

n t ’ 


or the ratio of energy that is not reflected 

a = 1 — r. 


(14) 


Here r is the statistical reflection co-efficient and is defined as 



(15) 


where II r is the power reflected. 

If there are other absorbing surfaces (such as seats) which are not part of the walls 
of the volume, the effects of the absorption of these other surfaces can be included in 
the modified room constant as 


# 0 / Non— Boundary Surfaces 

R' = S t a + £ S k a k . (16) 

fe=l 

If the volume is large, then the effect of absorption in the air must be added to the 
room constant. The air absorption term is 4 mV where m is the energy attenuation 
constant and V is the room volume. 
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Transmission Loss Theory and Validation 

3.1 Introduction 

Transmission loss ( TL ) is the classical measure of the installation-independent sound 
insulation properties of a panel, and it is defined in terms of the transmission co-efficient 

T ’ 1 

TL = 10 log (17) 

T 

In the modified room equation presented in Section 2, all walls or partitions are com- 
pletely characterized with transmission and reflection co-efficients, r and r. Much work 
has been done to predict and measure the transmission co-efficient. The reflection co- 
efficient has not received nearly as much attention, but it can be derived from the same 
equations that are used to predict the transmission co-efficient. 

The transmission loss model used in this study is based on an infinite panel theory. 
Although infinite panel theory does not take into account all the details of a fuselage 
sidewall or panel, it does predict the effects of several panel design parameters in 
a reasonable way [3]. The prediction of in situ transmission loss of a panel in its 
actual location in the fuselage or structure will require detailed theoretical modelling 
of the entire structure and installation. An infinite panel transmission loss theory, 
though simple in approach, will predict the effects of panel properties without detailed 
knowledge of the installation. 


3.2 Theory 

3.2.1 Definitions 

The transmission co-efficient r is 
power to the incident acoustic power 


Similarly, the reflection co-efficient r 
power to the incident acoustic power 


The absorption co-efficient is defined as the ratio of the acoustic power that is not 
reflected to the incident acoustic power or 

Cc — 1 — T. (20) 
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defined as the ratio of the transmitted acoustic 


T = 


Hi 

Hi' 


(18) 


is defined as the ratio of the reflected acoustic 


r = 


Ilr 

TT ' 


(19) 


I 

i 

i 

m 

t 

g 

| 
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3.2.2 Single Panel Transmission Loss 


The transmission loss equations for single and double wall infinite partitions are 
derived in many acoustics text books [3] [4], Since the TL equations are derived for 
a specific panel design, any deviation or added complexity in a panel would require a 
complete rederivation. The transmission loss equations for a single, infinite panel are 
given in this section. 

If all the equations for each type of medium are put in matrix form, the resulting 
equation is [ A]x = b where the variables associated with transmitted sound are in the 
unknown vector x, and the variables associated with incident sound are in the vector 
b. This allows the matrices for each portion of the panel to be cascaded together in 
order to get the complete panel matrix. This approach is identical to transmission line 
theory in electrical power transmission and is also called four-terminal (or four-pole) 
network theory. 

In the case of a single panel with a normally incident plane wave with geometry 
shown in Figure 1, the pressure and displacement variables can be written as 


= P ie ^- k ^ ( 21 ) 

Pr(x,t) = p r e^ t+k ^ (22) 

Pt(x,t) = P t e^- k ^ (23) 

f(M) = (24) 


where the p s are pressures, £ is displacement, k\ = ^ u> is the frequency in 

radians per second and the c’s are the speeds of sound in the media on both sides of 
the panel. The subscripts t, r and t denote incident, reflected and transmitted waves, 
respectively. The capital letters P and S are the pressure and displacement magnitudes 
of the acoustic wave at frequency u> and wave number k, respectively. The total pressure 
in medium 1 is the sum of the incident and reflected waves or 

Pi = Pi + Pr (25) 

and the total pressure in medium 2 is the pressure from the transmitted wave. At 
the interface between the air and the panel, there is continuity of displacement. The 
displacement continuity coupled with the linearized Euler’s equation 

du 

~VP = Po-Q t (26) 

(where p 0 is the ambient air density and u is the displacement of the acoustic medium) 
can be used to get the relationship between pressure and panel displacement for a 
periodic excitation, 

VP = PoU 2 £. (27) 
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Combining this relationship and equations (24) and (25) gives 


Pi~ Pr = JUPlCiE 

and 

Pt = JUP2C2E. 

The equation of motion of the panel is 

+ 77I+ < =Pi( 0 ,<) — p 2 ( 0 ,<) 


where m is the mass per unit area of the panel, 77 is the damping co-efficient of the 
panel and k is the stiffness of the panel. Combining equations (24), (25), (28), (29), 
and (30) produces 

AC 

jw5(jwm + 77-1 + pi ci + p 2 c 2 ) = ( z p + z h -| - z h )juE = 2 Pi (31) 

where the panel impedance is z p — jwm + 77 + and the impedance due to fluid loading 
on side i is Zf { = p{C{ . 

For plane waves, the transmission co-efficient is defined as 


mV 

2 pici 

and the reflection co-efficient is 


From equations (28), (29), (31), (32), and (33) r can be shown to be 


r — 4 


P ici \z h I 


and r can be shown to be 


P2C2 j z fl + Z p + Zf 3 


z h Zp + Zf 2 [ 

I Zfl + z P + z h\ 2 


As a special case consider a massive limp panel with the same air medium on each 
side or pc = p x Ci = p 2 c 2 m, k = 0 , and 77 = 0. In this case the transmission loss is 
given by 

r = C-^) 1 ( 36 ) 


which is the classical mass law formula for plane waves at normal incidence. 
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3.2.3 Matrix Formulation 


Consider the same panel again. At the interfaces between the panel and the acoustic 
medium, the variables of interest are the net pressure and the velocity or displacement. 
The velocities u on each side of the panel are the same or 


«i = u 2 (6 = 6)- 


(37) 


The equation of motion of the plate relates the pressures on either side of the panel 
and can be written as 

ZpU 2 (orl) = Pl- P2- (38) 

Equations (37) and (38) can be written in matrix form as 



(39) 


At this point equations (28), (29), and (39) can be used to derive r and r with the 
same results as before. If the panel design is modified by addition of another layer of 
different material for example, another equation could be written 


P 2 
u 2 


1 

0 



and by substituting equation (40) into (39) 



1 

0 

1 

0 



(40) 


(41) 

(42) 


In general the total panel matrix can be built up of several two-by-two matrices, one 
for each separate medium used to make the panel. These matrices are all multiplied 
together in the order of appearance in the panel. For example, a double wall panel 
with an air gap (Figure 2) would have an equation with three two by two matrices 
multiplied together such as 

[Total Panel Matrix^ 2 *^ = [ FirstPanel] { ' 2 * 2) [Air] ( ' 2x2) [SecondPanel] (2 * 2) (43) 


Once multiplied together the total panel matrix will be a non-symmetric complex 
two-by-two matrix and the algorithm for getting r and r will be identical for any panel 
design. This method of organizing the equations simplifies coding the equations for 
maximum flexibility in panel design. Also it is easy to include new types of materials 
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in the analysis because the only new equation that needs to be derived is the matrix 
equation which describes the behavior in the new medium. 

The transmission and reflection co-efficients can be calculated from their definitions 
(equations (32) and (33)), the total panel matrix and from Euler’s equation (26). The 
panel matrix equation may be written as follows: 



where the matrices [P«] and [P--\ are four-pole representation of panels 1 and 2 respec- 
tively, and the matrix [-A;,] defines the air gap between the two panels. 


The total panel equation can be rewritten in the following form: 

Tii T12 

T21 T22 

where p n and u n are the pressure and velocity at the panel surface in the last (nth) 
medium, respectively. The other equations lead to 




Pi = Pi + Pr (46) 

P;|*=o -Pr|*=o = Z h u 1 (47) 

Pn = Pt (48) 

Pt = z u u n . (49) 


After solving for the relationship between p t and pi and for the relationship between 
p r and pi, it is easy to get transmission and reflection co-efficients by substituting into 
their definitions, 


T 


Z fn 


( 


4 


\\Tn + % + z fl T 21 + 



(50) 


and 


T = 




(51) 


where Zf x is the acoustic impedance of the air at the surface of the panel on side 1, 
Zf n is the acoustic impedance of the air at the surface of the panel on the transmitting 
side of the panel (transmits into the nth medium) and is an element of total panel 
matrix. 
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The air impedances are pc for normal incidence, but for oblique incidence and other 
complicating factors such as air flow the impedance changes. The impedance of plane 
waves in air with flow at oblique incidence is [9] 


Z U = 


piCi 

sin 4>i( 1 + M cos/? cos 4>i) 


(52) 


where M is the Mach number of the flow, <f> is the angle of incidence with respect to a 
vector normal to the panel surface and (3 is the angle between the flow direction and 
the projection of the wave propagation on the panel surface. 

The matrix for an air gap between two panels is [5] 


cos (kidi cos(<j>i )) sin(fcjdi cos(^j)) 

j sin (kidi cos z fi cos (kidi cos (</>,)) 

where di is the depth of the air gap, Zj i is the acoustic impedance in the air gap, </>j is 
the angle of sound propagation with respect to a vector normal to the panel and k{ is 
the wave number of sound in the air gap. 

The matrix for a fiberglass layer is given by [5]: 


(53) 


cosh(qd m ) jZ mx sinh(qd m ) 
jsinh(qd m )/Z mx cosh (qd m ) 


(54) 


where d m is thickness of the fiberglass blanket. The propagation constant q and the 
characteristic acoustic impedance Z mi of the fiberglass material in the x-direction are 
given by [6]: 


9 = Trn 



jksintpi 

7m 


(55) 



(56) 


The propagation constant 7 m and the characteristic impedance Z m may be presented 
as a function of the dimensionless parameter 


7m — Ct rn -f- j/3 m 


where 


Zm — Rm + ]X m 
a„ = -[i 0.189(^)-° 595 ] 

C Jtii 

(3 m = -[l + 0.0978(^)-°' 700 ] 

C xti 


(57) 

(58) 

(59) 

(60) 
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( 61 ) 


R m = pc[i + o.om(^-y° 7M } 

ill 

= -pc[0M70(^-r° 732 } (62) 

ill 

where Ri is the flow resistivity (rayls/m) and p is the gas density in the ambient 
medium. 


Comparison with experimental data has shown that as a conservative approach, the 
range over which the above empirical equations for the semirigid materials are valid 
is [6,7]: 0.01 < < 0.1. In addition, the lower limiting frequency is determined by 

d m < A m , where A m is the wavelength of sound within the fiberglass blanket (A m — 




This method of setting up the transmission equations is very attractive due to the 
flexibility available in panel design and also to the ease of integrating new materials or 
more detailed analysis of old materials into the transmission loss prediction program. 


3.2.4 Mechanical Bridging Effects 


When a panel is made up of multiple layers and air gaps, there are usually frames 
holding the panels in place. This mechanical connection between the panels provides 
a flanking path for the energy to travel. Thus far only the airborne portion of the 
path has been considered. In this section a lumped parameter model of the structural 
path through the frames is developed. The frames are assumed to have periodic, 
parallel line connections on the panels. This type of connection is like the frames in 
a standard aircraft sidewall or like the studs in a wall in standard residential housing 
construction. This development is based on the four-pole (or four terminal) network 
theory [8]. The four-pole modelling approach permits a complete matrix formulation 
of a typical aircraft sidewall panel. The effects of fiberglass insulation can also be 
incorporated. 

For simplicity, a connected panel-airgap-panel system, as shown in Figure 3, is 
considered here. A structural path between the panels is provided by the frames. Two 
assumptions about the frames are: 


• the dynamic behavior of each frame is independent of any other frame 

• the frames translate as rigid bodies in their own planes. 

The first panel is excited by the acoustic wave driving a bending wave along the 
panel. When the wave hits a frame, the frame scatters the wave. Also the frame is 
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driven by the deflection of the first panel. The frame excites the second panel which 
radiates the sound. The structural sound path is assumed to be independent of the 
acoustic path. 

Neglecting the fluid loading effects, the bending wave velocity of the first panel is 

Mi ~ — , (63) 

jwm, 

where £1 is the displacement of the first panel, pi is the pressure incident on the panel, 
mi is the mass per unit area of the first panel and u; is the frequency in radians per 
second. The velocity of a frame (fine mass) when driven by a bending wave in a panel 
is 

> _ M-Zf, 

3 f (Z Fl + jwrrif + Zf 2 Y 
where is the displacement of the frame, m/ is the mass per unit length of the frame, 

Z Fi = 2Z) t 1/4 u; 1/2 mf /4 (l + ;) is the fine force impedance of panel i and D is the bending 
stiffness of panel i. 

The bridging impedance due to the mass of the frame and the line force impedances 
of the panels is therefore: 

Zb = Z Fl + jurrrif + Z Fl (65) 

In the four-pole representation of the connected panel system the bridging path 
(through the frame) can be considered as a shunt element, Zb, parallel to the airgap 
impedance, Z a . The force at the input point of the parallel connection is same as that 
at its output terminal. The total panel matrix for calculating transmission loss of such 
a panel configuration may now be written as follows: 



where the \(3\ matrix describes the four-pole parameters of the parallel connection. 
These parameters can be determined using matrix techniques [8]. 

The matrix equation for a general case of n impedances connected in parallel may 
be written as follows: 



where 



(67) 

( 68 ) 
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( 69 ) 


(70) 

(H) j 

| 

(72) | 

\ 

(73) | 

(74) | 

| 

3.3 Validation { 

The transmission loss prediction method presented in Section 3.2 was validated 
using laboratory data on test panels. The measured TL data for bare and treated j 

stiffened flat test panels (dimensions: 1.46m x 2.3m x 12mm), which are representative * 

of an aircraft fuselage sidewall, were used for the validation. j 

The predicted transmission loss of the bare test panel is compared with the labo- 
ratory test data in Figure 4. The TL of the test panel is underpredicted at low and 
mid-frequencies. Since the prediction model is based on the infinite panel theory, it 
does not account for the (two-dimensional) stiffening effect of frames and longerons on 
the test panel. 

The measured TL of the same test panel with fiberglass insulation (blanket thick- 
ness: 76.2mm) and a 1mm thick trim panel is compared with the predicted TL in 
Figure 5. The acoustic absorption effects due to fiberglass insulation were included in 

the predictions. The predicted TL is seen to be in good agreement with the test data. j 

1 

The comparisons of the predicted TL of a typical aircraft sidewall panel with labo- | 

ratory data therefore validates the TL prediction module. | 
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Figure 1: Geometry of normally incident plane wave on a single panel 
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Figure 3: Equivalent Impedance Representation of Panel-Airgap-Panel System with 
Structural Bridging 
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Figure 5: Predicted vs. Measured Transmission Loss Data of a Flat Sidewall Panel 
(.00127 m) with Fiberglass Insulation (0.0762 m) and Trim Panel (0.00076 m) 
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4 


Structure-Borne Noise Models 


4.1 Structure-Borne Noise Overview 

The structure-borne noise models were developed to predict the engine noise that 
is transmitted from the engine to the fuselage through a structural path (as opposed 
to an airborne path which impinges on the structure). The vibration of the engine 
is transmitted through the pylon (or wing) to the fuselage which radiates sound into 
the cabin. The inputs necessary to make this prediction are the forces of the engine 
acting at the engine mounts and the physical or statistical parameters which describe 
the pylon/fuselage region. 

The structure-borne path was separated into two regions. The first region (Region 
I) covers the path of the vibrational energy from the engine to the vibrational energy 
in the fuselage. The second region (Region II) covers the radiation of the fuselage 
vibrational energy into acoustical energy in the cabin. The energy radiated into the 
cabin is input into the energy balance equation as an additional energy source. 

For both regions, there are two alternative analysis procedures which may be used. 
In the first region the vibrational energy of the fuselage may be determined with either a 
simple SEA model of the pylon(or wing) /fuselage connection or a user defined transfer 
function. In the second region, the acoustical radiation from the fuselage may be 
predicted with either the radiation efficiency formula for flat plates or a user defined 
transfer function. 

4.2 SEA Model 

Since the Statistical Energy Analysis (SEA) approach is based on the overall statis- 
tical properties of energy flow, complex structures may be modelled by simple idealized 
structural elements with equivalent statistical properties. When using the SEA tech- 
nique, the dynamic system being studied is divided into a set of connected substruc- 
tures. SEA substructures are typically simple structural/fluid element representations, 
such as beams, plates, acoustic spaces etc., for sections of the modelled system. Each 
substructure is then divided into subsystems which are groups of modes having similar 
dynamic characteristics. The power balance equations are then set up by balancing the 
time average net vibratory power input to each subsystem with the time average power 
dissipated within the subsystem and the power transferred to other subsystem(s). Vi- 
bratory energy flow between subsystems is controlled by equations analogous to the 
equations for steady state heat or fluid flow. Energy flow between subsystems is pro- 
portional to the difference in subsystem modal energies. The input excitation for the 
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SEA model is applied as the power input to the specified subsystem. The power balance 
equations are then solved for the modal energies within the subsystems. The response 
variables of interest for the specified subsystem are determined from the modal energies. 

In order to keep the prediction scheme simple and general in scope so that it is 
appropriate for both wing and pylon mounted engines, the pylon/fuselage region is 
modelled as a two element SEA system. Each element of the SEA system is repre- 
sented by an equivalent flat plate. An equivalent flat plate is defined as a flat plate 
with statistical properties (such as modal density) the same as the pylon or fuselage 
structure. The pylon’s equivalent plate is rigidly connected perpendicularly to the fuse- 
lage’s equivalent plate. Engine vibration data are input to the pylon element and the 
SEA power flow equations are solved to give the space-averaged velocity in the fuselage 
which will be the input to the radiation model in Region II. 

The equations for power balance in SEA involve certain parameters defining the 
energy flow between subsystems, called “coupling loss factors”. The coupling loss 
factor is a measure of intermodal forces between the interacting multimodal systems, 
averaged over frequency and over the modes of the interacting systems. In addition, 
each SEA subsystem is statistically described by the following parameters: mode count 
iV(a>), wavenumber fc(ar), and damping loss factor rj. 

The steady state power balance equations for the case of two coupled subsystems 
are given by (see Figure 6) [11] 


nr = nf 9 + n 12 
nr = nf 9 - n 12 . 

The power dissipated in each subsystem is given by 

nf 9 = 


(75) 

(76) 

(77) 


and the net power flow from subsystem 1 to 2 is given by 

n 12 = ni 2 ' - n 21 ', (78) 

where H*/ is the power transmitted from subsystem i to subsystem j. nj n denotes the 
input power to subsystem 1 and E m represents the average subsystem modal energy, 


where E is the total energy of the subsystem. 

The net power transmitted from subsystem 1 to subsystem 2 is proportional to the 
difference in subsystem modal energies and may be expressed as 
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where E™ and E ™ are the average modal energies of subsystems 1 and 2 respectively. 

In matrix form, the power balance equations are 

{Vi + 7712 )^ 1 ( 0 ;) -V 2 iN 2 (v) 

(jJ 

-7712^(0;) (772 + 7721)^2 (w) _ 

There are essentially four parameters that enter into the SEA prediction scheme: 
(1) mode count, (2) coupling loss factor, (3) dissipative loss factor, and (4) input 
power. The mode count, N{oj ), of a subsystem is defined as the number of resonant 
modes within a given frequency band. The mode count for transverse vibrations of a 
homogeneous rectangular panel is given by [11] 



Ni(u>) = 


ojS, 


Pi 


4nK Pi Cu 


(82) 


where S Pi is the surface area of plate i, K Pi = is radius of gyration, c/ ; = 
is longitudinal wave speed, is the Poisson’s ratio for plate material and pi is the 
material density (mass per unit volume). 


The coupling loss factor describes the power flow from one subsystem to another. 
In the present case, the coupling is in the form of a mechanical connection between two 
structural elements. In the case of two plates connected at right angles, the coupling 
loss factor, 7712 , is related to the random incidence transmission co-efficient, t&, by [12] 


Vij = 


2^,- 

TtkiS Pi 


'ij > 


(83) 


where Lij is the length of the coupling fine and ki is the bending wavenumber in plate 
i. 


The random incidence transmission co-efficient, is [13] 


2 2.754X 

T ij w IX 5 / 4 + A" -5 / 4 ] 2 ' 1 + 3.24X 


(84) 


where X = h-i/hj , h, is the thickness of panel i and hj is the thickness of panel j. 

The predicted modal energy of a subsystem can be transformed into physical pa- 
rameters, such as spatially averaged vibration velocity, < v > by 
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( 85 ) 


. |2 E m N{u>) 
< v 2 > = v ' 


PsS p ’ 

where p, is the mass per unit area of the structural element. 


The mechanical input power injected into a structure by external sources is one 
of the parameters that is presumed known in the SEA prediction scheme. The input 
power at a given frequency or in a band of frequencies may be defined in terms of the 
input mobility function as follows 


n 


in 

1 


= ±\F\*Re{Y(w, I ,)} = \\V\* 


|i>, *<,)!’ ’ 


( 86 ) 


where Y (a>, x a ) is the input mobility of the structure at the driving point x a and may be 
defined as the complex ratio of the velocity vector V(w,x 0 ) to the force vector F(u,x 0 ) 




V{v,*o) 

F{u},x 0 y 


(87) 


and Re{..} indicates the real part of the argument function. 


The input mobility and the input force or velocity are used to calculate the input 
power as in equation (86). The physical parameters of the equivalent plate are used in 
equations (82)-(84) to calculate the terms in the SEA matrix equation (81). By solving 
equation (81), the power in the fuselage is found and may be used to calculate the 
space average velocity by equation (85). The space average velocity is then input to 
the radiation module (Region II). 


4.3 Radiation Efficiency Model 

The original work which defined the radiation efficiency of stiffened panels was done 
by Maidanik [14]. The form used in this analysis was presented by Ver and Holmer 
[3] . The acoustical power, II, radiated by a flat plate subjected to a random incident, 
uniform acoustical loading is given by 

n(a>) = e(u >)pcS < |v(o?)| 2 > (88) 

where e is the average radiation co-efficient, as a function of frequency, p is the density 
of the medium the panel is radiating into, c is the speed of sound in the same medium, 
S is the area of the radiating surface, and < v(u >) > is the space averaged velocity of 
the plate at frequency u>. 

The average radiation co-efficient, e, is defined differently in four separate frequency 
regions. These regions are: well below coincidence, below coincidence, at coincidence 
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and above coincidence. The coincidence frequency is: 





( 89 ) 


where D = * s the bending stiffness, c is the speed of sound in the radiating 

medium, p , — ph is the mass per unit area of the panel, p is the density of the panel, 
h is the thickness of the panel, E is Young’s modulus of the panel and v is Poisson’s 
ratio of the panel. 

In the four frequency regions, the average radiation co-efficient e is defined as follows: 


and 


where 


e(u>) = - w 2 when w < u> c and w < 


2s/YW 


7 r 2 c 2 


h^(MAX(l Xi l y )) 

2\/l2D 


p ^ 

(l \ 1/2 f l \ 1/2 

e(w) ~ I I M- when u) = w c ; 
e(u>) « ^1 when a? > w c ; 


9i(<r) = 


9i{v) ~ 


8 (l — ) * r /rv 

ir< a(l-<T 2 )‘/ 2 l f & < U> c /2 

0 if w > w c /2 ' 


[(l-^)ln(S) + 2^ 


4 7T 2 


(90) 

(91) 

(92) 

(93) 

(94) 

(95) 


(96) 

(97) 


(1 - tr 2 ) 3 / 2 

P = 2[/*(l + N ax ) + f y (l + iV, y )], 

A c = is the bending wavelength, N sx is the number of lengthwise stiffeners, N, y 

is the number of widthwise stiffeners, l x is the length of the panel and l y is the width 
of the panel. 


The radiation co-efficient e is calculated from the physical plate parameters of the 
radiating surface using equations (90) to (97). The radiation co-efficient together with 
the space-averaged velocity, equation (85) from Region I is used in equation (88) to 
calculate the power radiated from the plate. The radiated power is then input into the 
room equation as an interior power source. 
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4.4 Transfer Function Models 


The transfer function approach brings flexibility to the analysis of the structure- 
borne path. Experimental data or theoretical data from an alternative analysis may be 
included by the use of transfer functions. The transfer function equation for the first 
path (Region I) which is from the engine force to the fuselage vibration is 

F( u) x Hi{ u) = v(v), (98) 

where F is the engine force on the pylon, is the frequency in radians per second, Hj is 
the transfer function for Region I and v is the velocity of the fuselage. Both the engine 
force F and the transfer function Hj are input by the user and the fuselage velocity v 
is the output to be used in Region II. 

The transfer function equation for the second path (Region II) which is from the 
the fuselage vibrational velocity to the radiated acoustical power is 

v(uj) x Hn(u>) = II(a>), (99) 

where v is the velocity of the fuselage, w is the frequency in radians per second, Hu is 
the transfer function for Region II and II(u;) is the acoustical power radiated into the 
cabin. In this region, the fuselage velocity v is the result from Region I, or from SEA 
model, equation (85) the transfer function Hn is input by the user, and the acoustical 
power n(o>) computed by equation (99) is input into the room equation as an alternate 
source of energy. 
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Figure 6: Power Flow Between Two Sets of Linearly Coupled Subsystems 
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5 Boundary Layer Noise Model and Validation 

5.1 Prediction Model 


The pressure fluctuations in the turbulent boundary layer caused by the airflow past 
the aircraft are one of the main sources of cabin noise. The turbulent boundary layer 
prediction scheme used in ATDAC is based on empirical formulae derived by Lowson 
[15]. 

The overall root mean square pressure fluctuation of the spectrum, p rmi , is given 
by the empirical equation 

Prm, _ 0-006 

q ~ 1 + 0.14 M 2 , 

where q = \p 0 U 2 is the dynamic pressure head, p 0 is the air density, U = cM is the 
aircraft velocity, c is the speed of sound and M is aircraft Mach number. The power 
spectral density is defined by another empirical equation 



p(w ) 2 = 


PL, 

u; 0 [l + (a;/u; 0 ) 2 ] 3 / 2 


( 101 ) 


where w is the frequency (rad/sec) and u> 0 = 8U/8 is a characteristic frequency for nor- 
malization purposes. The boundary layer thickness 8 is given by an empirical formula 
suggested by Bies [16] 


- = 0 . 37tf ei - 0 - 2 

X 



Re^ 

6.9 x 10 7 



(102) 


where x is the distance from the nose of the aircraft, Re x = XJxjv is the Reynolds 
number and v is the kinematic viscosity. The one-third octave pressure spectrum is 
calculated from the power spectral density by integrating over each one-third octave 
band. The power spectral density is assumed to be smooth and fairly slowly varying 
for integration over each one-third octave band. 


5.2 Validation 

This empirical model was compared with two sets of experimental data. All the data 
were non-dimensionalized to simplify the comparison. The magnitud e of the meas ured 
boundary layer pressure was non-dimensionalized by dividing by q^j (T ^o.i4M 2 ) an< ^ 
the frequency was non-dimensionalized by multiplying by The empirical prediction 
curve collapses to a single curve when non-dimensionalized in this way. The non- 
dimensionalized data in Figure 7 for the MD-UHB Demonstrator (which was limited 
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to frequencies below 2000 Hz before normalization) display a lot of scatter. The data 
in Figure 8 for the Boeing 737 [17] aircraft show similar results but extend to higher 
frequencies. In the mid to low frequency region the empirical model seems to provide 
an upper bound to the range of measurements, but in the higher frequency region 
the empirical model drops off faster than the measured data. The empirical model is 
conservative in the mid to low frequency regions and for many noise problems these 
are the important frequencies. 
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Figure 7: Comparison of Predicted (— ) and Measured (o, •) Turbulent Boundary Layer 
Noise. Measured Data from MD-UHB Demonstrator Aircraft. 
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Figure 8: Comparison of Predicted (-) and Measured («— , f , J., •) Turbulent Boundary 
Layer Noise. Measured Data from a Boeing 737 Aircraft. 
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6 Cabin Noise Prediction Validation 


The validation of ATDAC involved comparison of predicted transmission loss, bound- 
ary layer noise and cabin noise levels with measured data. The transmission loss and 
boundary layer noise prediction modules were validated with experimental data in 
Sections 3.3 and 5.2 respectively. In this section, the overall ATDAC interior noise 
predictions are validated using flight test data for the MD-UHB Demonstrator aircraft. 
The Demonstrator aircraft is a modified MD-80 series aircraft. The MD-80 aircraft, 
with an overall length of 147.8 ft and a wing span of 107.8 ft, is designed to carry 155 
passengers. The MD-80 aircraft is normally powered by two Pratt & Whitney JT8D- 
209 engines. The aircraft (ship 909) was modified by replacing the left engine with a 
General Electric proof-of- concept Unducted fan (UDF) engine. 

The validation of ATDAC as an aircraft cabin noise prediction method required 
complete modelling of the MD-UHB Demonstrator aircraft. The details of aircraft 
cabin, noise sources and material modelling are given in the User’s manual [1]. The 
aircraft cabin was divided into 47 volumes and 111 exterior surfaces. There were 74 
interior partitions, such as the floor and pressure bulkhead, which separated different 
volumes. In the case of two adjoining volumes with no material surface joining them, 
an air partition was used. The geometrical details, such as surface areas, global loca- 
tion of each volume and partitions joining different volumes, were generated using the 
geometry module and stored in a geometry file. The transmission loss and reflection 
co-efficients of the fuselage surfaces and interior partitions were estimated using the 
transmission loss module and were stored in a panel database file. The acoustic ab- 
sorption characteristics of the interior cabin surfaces were contained in absorption files. 
A material file assigned the acoustical characteristics of the surfaces and partitions 
according to the geometrical description of the aircraft cabin. 

The source file module prepared an output source file which contained information 
about all the noise sources including exterior boundary layer noise, airborne noise due 
to the propfan and structure-borne engine noise. The UHB engine used for these flight 
tests was a General Electric UDF engine with two rows of highly swept, counter rotating 
propeller blades. The flight test data from an 8x8 engine configuration (consisting of 
eight blades on both the forward and aft rotors) was used for the ATDAC validation 
in this report. Interior noise spectra (see Figure 9) in the aft cabin, measured in flight, 
were found to contain several tones, superimposed on broadband noise. These tones 
correspond to the blade passage frequencies of the propeller and other engine sources. 
Tones at the blade passage frequency (BPF) and harmonics are generated by each 
propeller. These propeller harmonics are labelled as 2BPF, 3BPF etc. An additional 
tone labeled UNI is associated with the rotational speed of the low pressure stage of the 
General Electric UDF engine. The broadband noise is mainly due to exterior boundary 
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layer noise. 

A knowledge of the main noise transmission paths into the aircraft cabin is very 
useful in modelling various noise sources. Extensive laboratory and flight test inves- 
tigations conducted on the MD-UHB Demonstrator aircraft provided measured data 
for ATDAC prediction model and information for modelling noise transmission paths. 
The following noise transmission paths were modelled: 

• an airborne noise path through the cabin sidewall for propeller tones and broad- 
band boundary layer noise component 

• a structure-borne transmission path through the engine pylon and the aft fuselage 
for propeller and engine tones 


• a combined airborne/structure-borne path for propeller tones: the tail section 
is set into vibrations due to acoustic energy from the propeller blades; the aft 
pressure bulkhead, in turn, is excited by the tail surface and radiates energy into 
the cabin. 

Measured flight test data for the MD-UHB Demonstrator aircraft provided input 
levels for airborne and structure-borne noise sources [18]. The main structure-borne 
noise source for cabin noise prediction was simulated using the measured vibration 
levels at the left aft mount of the UHB engine. Figure 10 shows measured acceleration 
levels at aft engine mount location for BPF and UNI tones. 

The structure-borne sound transmission from the engine pylon was modelled using 
the SEA approach described in Section 4. The input power due to structure-borne noise 
sources was estimated using the drive-point input mobility formulation. In addition to 
boundary layer noise, propeller tones at blade passage frequencies were simulated as ex- 
terior airborne noise sources for the aft fuselage. Figure 11 shows exterior noise spectra 
on the aft fuselage due to the 8x8 propfan engine. A combined airborne/structure- 
borne path through the aft fuselage (tail surface) and aft pressure bulkhead was simu- 
lated using the transfer function approach. Figure 12 shows BPF vibration and noise 
levels as a function of rotor speed for a fuselage accelerometer in the forward prop plane 
and an interior microphone in the unpressurized section aft of the bulkhead. 

The cabin noise levels due to exterior airborne noise source, i.e. propeller tones 
1BPF and 2BPF, and structure-borne noise sources, i.e., engine tone UNI and pro- 
peller tone 1BPF, were predicted separately. The effect of airborne noise source due to 
exterior boundary layer noise was not included. The predicted airborne and structure- 
borne cabin noise components are compared with measured flight test cabin noise levels 
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in Figure 13. The predicted cabin noise levels (spatially-averaged levels in Row 6 seats) 
shown in Figure 13 compare well with the measured data. It was found that the com- 
bined airborne/structure-borne path, e.g., via aft fuselage (tail section) and pressure 
bulkhead for 1BPF propeller tone (at 168 Hz) was more dominant than the airborne 
path. 

A comparison of ATDAC predicted cabin noise levels with flight test data for aft 
cabin locations (Row 6 seats) on the Demonstrator aircraft is shown in Figure 14. 
The predicted interior noise levels are in good agreement with those measured during 
flight. The cabin noise levels above 200 Hz 1/3-octave band are mainly controlled by 
the turbulent boundary layer noise. The structure-borne cabin noise component due 
to UNI engine tone in the 200 Hz 1/3-octave band was overshadowed by the airborne 
contribution from exterior turbulent boundary layer noise. The airborne contribution 
from 2BPF propeller tone in the 315 Hz 1/3-octave band was similarly masked by 
the boundary layer noise component. The overestimation of cabin noise levels in the 
frequency range of 200 and 800 Hz may be attributed to the following: (1) higher 
exterior boundary layer noise spectrum levels and, (2) lower transmission loss of the 
fuselage. The predicted cabin noise levels are in very good agreement with the measured 
data at both low and high frequencies. 
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Figure 10: Aft Engine Mount Acceleration Levels - 8x8 UHB Engine 
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Figure 13: Comparison of ATDAC Predicted Airborne and Structure-Borne Cabin 
Noise Levels with Flight Test Data for MD-UHB Demonstrator Aircraft. (Altitude: 
35000 ft, 0.76 Mach, 1260 RPM) 
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Figure 14: Comparison of ATDAC Predicted Cabin Noise Levels with Flight Test Data 
for MD-UHB Demonstrator Aircraft. (Altitude: 35000 ft, 0.76 Mach, 1260 RPM) 
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7 Conclusions 


An analytical model for prediction of noise levels inside an aircraft cabin based on 
power balance, the acoustic room equation and transmission loss concepts has been 
developed. The theoretical basis for the computer program ‘Acoustical Theory for 
Design of Aircraft Cabins (ATDAC)’ is presented in this report. This program is a 
set of modules based on validated algorithms for prediction of turbulent boundary ! 

layer noise, transmission and reflection co- efficients of aircraft-type panels, structure- j 

borne noise and interior noise levels. In ATDAC, these modules form a system to 
predict the aircraft interior noise levels from jnformation about the acoustic sources, 
aircraft geometry and design. The details of the implementation of these algorithms are 
described in the following Report: “Interior Noise Prediction Methodology- ATDAC 
User’s Manual.” [1] j 

Flight and ground test interior noise and vibration data were used to validate the ? 

code’s predictions. The predicted transmission loss for selected panel configurations j 

were compared with laboratory test data. The boundary layer noise prediction model | 

was validated using flight test data. The ATDAC predicted cabin noise components J 

due to dominant airborne and structure-borne paths were compared with the measured \ 

data. The ATDAC interior noise prediction scheme was validated using flight test data ] 

for the MD-UHB Demonstrator aircraft. The ATDAC predicted cabin noise levels are ; 

in good agreement with the flight test data. ? 

From the results presented in this report, it may be concluded that ATDAC provides 
a useful interior noise prediction methodology for aircraft cabins. A noise control I 

designer can use the ATDAC program to compare the effects of different noise control 5 

treatments and thereby develop an aircraft configuration which meets interior noise i 

criteria. 1 
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with advanced turboprop (ATP) aircraft where there is a tonal, low frequency noise problem. 

Prediction of interior noise levels requires knowledge of the energy sources, the transmission paths, and the relationship . 
between the energy variable and the sound pressure level. The energy sources include engine noise, both airborne and 
structure-borne; turbulent boundary layer noise; and interior noise sources such as air conditioner noise and auxiliary power, 
unit noise. Since propeller and engine noise prediction programs are widely available, they are not included in ATDAC. 
Airborne engine noise from any prediction or measurement may be input to this program. This report describes the theory 
and equations implemented in the ATDAC program. 1 
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